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Analysis of hadron production experimental data allows to understand the prop- 
erties of the dense matter fireball produced in relativistic heavy ion collisions. We 
interpret the analysis results and argue that color deconfined state has been formed 
at highest CERN-SPS energies and at BNL-RHIC. 



1 Quark-Gluon Plasma Phase Boundary 

The signature of quark-gluon plasma is hard to identify, since the final state 
observed in the laboratory always consists of the same particles, irrespective of 
the transitional presence of the deconfined state. What changes is the detailed 
composition_jOf the observed produced particle abundance. We consider here 
strangeness El, and the pattern of strange hadron production at SPS and RHIC 
which offer, as we shall see, a compelling evidence for quark-gluon plasma 
formation. 

An expanding fireball of quark and gluons breaks up into final state 
hadrons in conditions _which differ from the equilibrium transformation ex- 
plored in lattice QCD □. To understand this behavior we need to understand 
the pressure of quark-gluon plasma: 

We have inserted here the appropriate quark (with ?it = 2) and gluon degen- 
eracy. The interactions between quarks and gluons manifest their presence 
aside of the vacuum structure effect B, in the three coefficients q ^ 1, see H: 

Cl = 1 - — ^ + • • • , C2 = 1 - — ^ • • • , C3 = 1 ^ + • • • . (2) 

Ait ^ ZItt ^ . . 

While higher order terms in the perturbative expansion have been obtained, 

they suggest lack of convergence of the expansion scheme formulated today. 

To reproduce the lattice results, near to T = Tc, the bag constant, B = 

0.19GeV/fm^ is used. 



Q^^+^^is^^'^^) +15^ 
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The condition P = Pp — B ^ 0, including the effect of motion reads 

S = Pp + (Pp + ep)^, .= ^^, (3) 

where subscript p refers to particle pressure and energy, and where we intro- 
duced the geometric factor k which characterizes the angular relation between 
the surface normal vector and flow direction. The motion of quarks and giu- 
ons exercises additional pressure and thus the QGP phase remains the matter 
phase even below the static equilibrium transition temperature. 

A fireball surface region which reaches condition Eq. (H) and continues 
to flow outwards must be torn apart. This is a collective instability and 
the ensuing disintegration of the fireball matter should be very rapid. A 
rapidly evolving fireball which supercools is in general highly unstable, and we 
expect that a sudden transformation (hadronization) into confined matter can 
ensue in such a condition. The situation we described could only arise since 
the vacuum pressure term is not subject to flow and always keeps the same 
value. An exploding QGP fireball only contained by the vacuum supercools 
to T ~ 145 MeV. Hadrons are produced well below the equilibrium transition 
temperature T ~ 160-170 MeV. 

2 Statistical Hadronization 

An extra reaction step 'hadronization' is required to connect the properties 
of the deconfined quark-gluon matter fireball, and the experimental appara- 
tus. In this process, the quark and gluon content of the fireball is transferred 
into ultimately free flowing hadronic particles. In hadronization, gluons frag- 
ment into quarks, and quarks coalesce into hadrons. Color 'freezes', quark- 
giuon plasma excess entropy has to find a way to get away, so any additional 
production is hindered. It is far from obvious that hadron phase space (so 
called 'hadronic gas') be used consistently to describe the physics of thermal 
hadronization, and we establish now the consistency criterion assuming that 
entropy production is small, or even null, in hadronization of entropy rich 
thermal QGP into entropy poor hadron phases space. 

Using the Gibbs-Duham relation for a unit volume, e + P = Ta + fii' , 
the instability condition of dynamical expansion, Eq. (^, takes the form: 

1 — f^h V (Th/V 1 — Wc/ 

Using extensive variables we obtain 
E 



S 



h 



For RHIC, we have (5r|h < 0.4 MeV, considering that < 40 MeV and 
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S/b > 100; at top SPS energy, we have ^ 200-250 MeV and S/b ~ 25-45, 
and thus, 5T\^~b MeV. 

The usefulness of Eq. (||) comes from the observation that it impHes: 

f >TW- (6) 

This is a near equahty since the geometric emissivity factor k is positive and 
small, especially so at RHIC, and < 1/3. 

The Gibbs-Duham relation also imphes: E/S + PV/S = T + 6T > T . 
One would think that the PV term is small, since the pressure is small as 
the lattice calculations suggest, e/P — > 7. However, this in principle can be 
compensated by large volume of hadronization. Since the volume is a directly 
measured quantity, this hadronization requirement can not remain unnoticed 
for long. In fact, the HBT results today place a very severe constraint on the 
emitter size of the pion source. With realistic volume size, the PV term is 
neghgible. In this case, we obtain just as in Eq. (||), E/S "^T . 

In Eq. (^), we have thus a necessary statistical hadronization constraint. 
In our experience, this is practically impossible to accomplish in chemical 
equilibrium statistical hadronization models. The reason is that, in such an 
approach, a rather high temperature T ~ 175 MeV is required to accom- 
modate the high intrinsic entropy content of hadron source, but this does 
not drive up sufficiently the energy content, since for pions in equilibrium 
E/S < T. What happens in such an approach is that a very large hadroniza- 
tion volume is required, in disagreement with the HBT results. We conclude 
that the combination of HBT result with fundamental properties of statistical 
physics requires chemical non-equilibrium statistical hadronization. 



3 Hadron Phase Space 

We GOiiisider a generalization of Fermi's statistical model of hadron produc- 
tion oB, and consider the yield of hadrons to be solely dictated by the study 
of the magnitude of the phase space available. 

The relative number of final state hadronic particles freezing out from, 
e.g., a thermal quark-gluon source, is obtained noting that the fugacity fi 
of the i-th emitted composite hadronic particle containing /c-components is 
derived from fugacities and phase space occupancies 7^: 

N,o,e-^-^^^.n^e-^-^^^l[j,X,. (7) 

In most cases, we study chemical properties of light quarks u, d jointly, 
though on occasion, we will introduce the isospin asymmetry. As seen in 
Eq. (M), we study particle production in terms of five statistical parameters 
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T, Aq, As, 7q, 7s- In addition, to describe the shape of spectra, one needs mat- 
ter flow velocity parameters, these become irrelevant when only total particle 
abundances are studied, obtained integrating all of phase space, or equiva- 
lently in presence of strong longitudinal flow, when we are looking at a yield 
per unit of rapidity. 

The difference between the two types of chemical parameters, Xi and 7^, is 
that the phase space occupancy factor 7^ regulates the number of pairs of flavor 
'z', and hence applies in the same manner to particles and antiparticles, while 
fugacity Xi applies only to particles, while A^^ is the antiparticle fugacity. 

The resulting yields of final state hadronic particles are most conveniently 
characterized taking the Laplace transform of the accessible phase space. This 
approach generates a function which, in its mathematical properties, is iden- 
tical to the partition function: 



C 



-Ei/T{ 



Y\_lkXk 



oclnZ«°. (8) 

Eq. (|^) does not require formation of a phase comprising a gas of hadrons, but 
is not inconsistent with such a step in evolution of the matter; it describes 
not a partition function, but just a look-alike object arising from the Laplace 
transform of the accessible phase space. The final particle abundances, mea- 
sured in an experiment, are obtained after all unstable hadronic resonances 
'j' are allowed to disintegrate, contributing to the yields of stable hadrons. 

The unnormalized particle multiplicities are obtained differentiating 
Eq. (H) with respect to particle fugacity. The relative particle yields are simply 
given by ratios of corresponding chemical factors, weighted with the size of the 
momentum phase space accepted by the experiment. The ratios of strange 
antibaryons to strange baryons of same particle type are, in our approach, 
simple functions of the quark fugacities. 



4 Strange Hadrons at SPS 

We expect, in sudden hadronization, chemical non-equilibrium at hadron 
freeze-out. Full chemical non-equilibrium description of particle yields is re- 
quired to arrive at a statistically significant description of the data with small 
error: 

(Ai?Lp)2 

It is common to normalize the total error by ths difference between the 
number of data points and parameters used, the so called 'dof (degrees of 
freedom) quantity. For systems we study, with a few degrees of freedom 
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Table 1. WA97 (top) and NA49 (bottom) Pb-Pb 158^ GeV collision hadron ratios 
compared with phase space fits. 



Ratios 


Ref. 


Exp. Data 






S/A 








0.099 ± 0.008 


0.096 


0.095 


H/A 








0.203 ± 0.024 


0.197 


0.199 


A/ A 








0.124 ± 0.013 


0.123 


0.122 










0.255 ± 0.025 


0.251 


0.255 


K+/K- 




i 




1.80± 0.10 


1.746 


1.771 


K-/7r- 






0.082±0.012 


0.082 


0.080 


K"Jb 




10 




0.183 ± 0.027 


0.192 


0.195 


h-/b 




n 




1.97 ± 0.1 


1.786 


1.818 


</./K- 




12 




0.145 ± 0.024 


0.164 


0.163 


A/p 


y 









0.565 


0.568 


p/tt 


all 


y 




0.017 


0.016 








1.6 


1.15 




N;p;r 




9;4;1 


9;5;1 



(typically 5~15), a statistically significant fit requires that x^/dof < 1. For 
just a few 'dof, the error should be as small as x^/dof< 0.5. The usual 
requirement ^ 1 is only applying for very large 'dof. 

Turning to the Pb-Pb system at 158 A GeV collision energy, we consider 
particle listed in table |l|, top section from the experiment WA97, for p± > 
0.7 GeV, within a narrow Ay = 0.5 central rapidity window. Further below 
are shown results from the large acceptance experiment NA49, extrapolated 
by the collaboration to full 4tt phase space coverage. The total error for 
the two result columns is shown at the bottom of this table along with the 
number of data points W, parameters 'p' used, and number of (algebraic) 
redundancies 'r' connecting the experimental results. For r 7^ 0, it is more 
appropriate to quote the total x^, since the statistical relevance condition is 
more difficult to establish given the constraints, but since x^/ {N—p—r) < 0.5, 
we are certain to have a valid description of hadron multiplicities. 

In second last column, the superscript 's' means that Ag is fixed by 
strangeness balance and, superscript '7q', in two last columns, means that 
7q = 7q = gm„/2Tf^ gjjed to maxlmlzc the entropy content in the hadronic 
phase space. The fits presented are obtained with latest_NA49 experimental 
results, i.e., have updated h~/6, newly published 4> yield Eil, and we predict the 
A/p ratio, b is here the number of baryon participants, and h^ ~ +K~ +p 
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Table 2. Upper section: statistical model parameters which best describe the exper- 
imental results for Pb-Pb data seen in Fig. 0. Bottom section: energy per entropy, 
antistrangeness, net strangeness of the full hadron phase space characterized by 
these statistical parameters. In column two, we fix As by requirement of strangeness 
conservation, and in this and next column we fix 7q = 7q. Superscript * indicates 
values which are result of a constraint. 









T [MeV] 


151 ± 3 


147.7 ± 5.6 




0.55 ± 0.05 


0.52 ± 0.29 


K 


1.617 ± 0.028 


1.624 ± 0.029 


As 


1.10* 


1.094 ± 0.02 


7q 




7^* = e"'/2'^'=1.6 


7s /7q 


1.00 ± 0.06 


1.00 ± 0.06 


E/b[GeY] 


4.0 


4.1 


s/b 


0.70 ± 0.05 


0.71 ± 0.05 


£;/S'[MeV] 


163 ± 1 


160 ± 1 


[s~.s)/b 


0* 


0.04 ± 0.05 



is the yield of stable negative hadrons which includes pions, kaons and an- 
tiprotons. We see, comparing the two columns, that strangeness conservation 
(enforced in second last column) is consistent with the experimental data 
shown, enforcing it does not change much the results for particle multiplici- 
ties. 

The six parameters (T, Wc, Aq, As, 7q, 7s) describing the particle abun- 
dances are shown in the top section of table |[ Since the results of the WA97 
experiment are not covering the full phase space, there is a reasonably precise 
value found for one velocity parameter, taken to be the spherical surface flow 
velocity Vc of the fireball hadron source. 

A value A^*^ ~ 1.1 characteristic for a source of freely movable straiige 
quarks with balancing strangeness in presence of strong Coulomb potential t3, 
i.e., with As = 1, is obtained. Since all chemical non equilibrium studies of the 
Pb-Pb system converge to the case of maximum entropy, we have presented 
the results with fixed 7q = 7q = e™x/2Tt_ ^pj^g large values of 7q > 1 confirm 
the need to hadronize the excess entropy of the QGP possibly formed. This 
value is derived from both the specific negative hadron hr /b abundance and 
from the relative strange hadron yields. 

The fits shown satisfy comfortably the statistical hadronization constraint 
that E/S > T discussed in section ||. We see also that near strangeness 
balance is obtained as result of the fit. 
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One of the interesting quantitative results of this analysis is shown in the 
bottom section of table ||: the yield of strangeness per baryon, s/6 ~ 0.7 . As 
we will show in Fig.|l], the expected equilibrium yield is s/5|*^*^^ ~ 1.4, and 
thus we have 7^^^ — 0.5 at SPS. Since the occupancy values in table || are 
derived from hadron phase space, we thus find 7^*^/7?^^ — 3. 

5 Strangeness at RHIC 

In the likely event that the QGP formed at RHIC evolves towards strangeness 
chemical equilibrium abundance, or possibly even exceeds it, we should expect 
a noticeable over occupancy of strangeness as measured in terms of chemical 
equilibrium final state hadron abundance. Because much of the strangeness 
is in the baryonic degrees of freedom, the kaon to pion ratio should appear 
suppressed, compared to SPS results. An even more penetrating effect of the 
hadronization of strangeness rich QGP at RHIC is the abundant formation of 
strange baryons and antibaryons Eil. This high phase space occupancy is one 
of the requirements for the enhancement of multistrange (anti)baxyon produc- 
tion, which is an important hadronic signal of QGP phenomena In partic- 
ular, we hope that hadrons produced in phase space with a small probability, 
such as 51, n, will be observed with a yield above statistical hadronization 
expectations, continuing the trend seen at SPS. 

The RHIC data we consider are were obtained in ^snn = 130 GeV run at 
the central rapidity region where, due to approximate longitudinal scaling, the 
effects of flow cancel and we can evaluate the full phase space yields in order 
to obtain particle ratios. We do not fit trivial results such as tt^/tt^ = 1, since 
the large hadron yield combined with the flow of baryon isospin asymmetry 
towards the fragmentation rapidity region assures us that this result will occur 
to a great precision. We also do not use the results for K* , K* since these yields 
depend on the degree of rescattering of resonance decay products. The data 
we consider has been reported by the STAR collaboration of Summer 2001, 
and where available is combined with data^PHENIX, BRAHMS, PHOBOS, 
for more discussion of the data origin, see E£l. We assume, in our fit in table ^, 
that the multistrange weak interaction cascading S — > A, in the STAR result 
we consider, is cut by vertex discrimination and thus we use these yields 
without a further correction. 

Wc first consider what experimental hadron yield results shown in table ^ 
imply about total strangeness yield in the RHIC-130 fireball. We begin with 
the yield of strange quarks contained in hyperons. We have, in singly strange 
hyperons, 1.5 times the yield observed in A, since remain unobserved. 
Also, accounting for the doubly strange which are half of the all S, and 
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contain two strange quarks, we have: 

= 1-5 • 0.059 + 2 • 2 • 0.195 • 0.059 = 0.133 . 

n 

Allowing for the unobserved Q at the theoretical rate, this number increases 
to (s)Y/h~ = 0.14. Repeating the same argument for antihyperons the result 
is 0.10. s and s content in kaons is four times that in Ks and thus the total 
strangeness yield is 

^=0.76, 
h 

with 32% of this yield contained in hyperons and antihyperons. This amounts 
to about 1.5 fold enhancement compared to highest SPS energies, with 
(multi)strange hyperons and antihyperons being a very important component. 

We now discuss the fits of RHIC-130 results. In the last column in table ^, 
the chemical equilibrium fit, the large originates in the inability to account 
for multistrange S, S. Similar results are presented in Ref. which work 
does not include multistrange hadrons. This equilibrium fit yields E/S = 
159 MeV < T = 183 MeV contradicting the conditions we discussed in depth 
in section ^. On the other hand, the chemical nonequilibrium fits come out 
to be in near perfect agreement with data, and consistent with the QGP 
statistical hadronization picture: E/S ~ 163 > T ~ 158 MeV and 7s, 7q > 1 . 
The value of the hadronization temperature T = 158 MeV is below the central 
expected equilibrium phase transition temperature, and this hadronization 
temperatures at RHIC is consistent with sudden breakup of a supercooled 
QGP fireball. The inclusion of the yields of multistrange antibaryons in the 
RHIC data analysis, along with allowance for chemical non-equilibrium (7 ^ 
1), allows to discriminate the different (chemical equilibrium/nonequilibrium) 
reaction scenarios. 

The value of the thermal energy content E/b = 25 GeV, seen in table |[ is 
in very good agreement with expectations once we allow for the kinetic energy 
content associated with longitudinal and transverse motion. The energy of 
each particle is 'boosted' with the factor 7^ coshy||. For uj_ = c/a/3, we have 
7^ = 1.22. The longitudinal flow range is about ±2.3 rapidity units, according 
to PHOBOS results. To obtain the energy increase due to longitudinal flow, 
we have to multiply by the average, / dy\\ coshyn/yy sinh(2.3)/2.3 = 2.15, 
for a total average increase in energy by factor 2.62, which takes the full 
energy content to E'"/^ ~ 65 GeV as expected. This consistency reassures 
that we have a physically relevant fit of the data, which respects overall energy 
conservation. 

We analyze next the strangeness content, s/6 = 6, seen table ^. The fully 
equilibrated QGP phase space would have yielded 8.6 strange quark pairs per 
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Table 3. Fits of central rapidity hadron ratios for RHIC -^snn = 130 GeV run. 
Top section: experimental results, followed by chemical parameters, physical prop- 
erty of the phase space, and the fit error. Columns: data, full non-equilibrium fit, 
nonequilibrium fit constrained by strangeness conservation and supersaturation of 
pion phase space, and in the last column, equilibrium fit constrained by strangeness 
conservation, upper index * indicates quantities fixed by these considerations. 





Data 


Fit 


Fit 

s-s = 


Fit°q 

s - s = 


p/p 


0.64 ±0.07 


0.637 


0.640 


0.587 


p/h- 




0.068 


0.068 


0.075 


A/A 


0.77 ±0.07 


0.719 


0.718 


0.679 


A/h- 


0.059±0.001 


0.059 


0.059 


0.059 


A/h- 


0.042±0.001 


0.042 


0.042 


0.040 


S/S 


0.83 ±0.08 


0.817 


0.813 


0.790 


S-/A 


0.195±0.015 


0.176 


0.176 


0.130 


5-/A 


0.210±0.015 


0.200 


0.200 


0.152 


K-/K+ 


0.88 ±0.05 


0.896 


0.900 


0.891 


K-/7r- 


0.149±0.020 


0.152 


0.152 


0.145 


Ks/h- 


0.130±0.001 


0.130 


0.130 


0.124 


0/s- 




0.222 


0.223 


0.208 






0.257 


0.256 


0.247 






0.943 


0.934 


0.935 


T 




158± 1 


158± 1 


183± 1 


7q 




1.55±0.01 


1.58±0.08 


1* 


K 




1.082±0.010 


1.081±0.006 


1.097±0.006 


7s 




2.09±0.03 


2.1±0.1 


1* 


As 




1.0097±0.015 


1.0114* 


1.011* 


£;/6[GeV] 




24.6 


24.7 


21 


s/b 




6.1 


6.2 


4.2 


S/b 




151 


152 


131 


£;/S'[MeV] 




163 


163 


159 


xVdof 




2.95/(10-5) 


2.96/(10-4) 


73/(10-2) 



baryon at Aq = 1.08, as is shown in Fig.g below. Thus 7^^^ = 6/8.6 ~ 0.7, 
which is greater than the value 0.5 we found for the SPS energy range. Using 
the fitted value 7^ ^ = 2.1, we find again, just like we found in case of SPS, 
7^*^/7^*^^ — 3. The fact that the strangeness phase space in QGP is not 
fully saturated is, on a second carefuLlook, in qualitative agreement with 
kinetic strangeness theory predictions Ej, adjusting our study to the observed 
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RHIC-130 run conditions. 



6 Strangeness as Signature of Deconfinement 

We consider the ratio of equilibrium strangeness density, arising in the Boltz- 
mann gas hmit, to the baryon density in a QGP firebaU. To first approxima- 
tion, perturbative thermal QCD corrections cancel in the ratio. For nis — 200 
MeV and T = 150 MeV, we have: 

b~^' lnAq + (lnAq)3/7r2 ^ 
The relative yield s/b is mainly dependent on the value of Aq, with only a 
slow temperature dependence contained in the coefficient 0.7 in Eq. (|o|). The 
light quark fugacity Aq is usually independent of the strategy of data analysis, 
as is also seen in table |^. 

At top SPS energy where Aq ~ 1.6, we see in Fig. |] that the equilibrium 
strangeness yield is at 1.5 strange pairs per participating baryon. The actual 
experimental yield, 0.7, see table ^, is half as large as in an equilibrated 
QGP, but is 2.5 times the yield in p-p reactions. At the RHIC 130 GeV run, 
the value Aq — 1.08, see table |3| and the specific strangeness yield in a QGP 
fireball at equilibrium is an order of magnitude greater than currently observed 
at SPS top energy. As discussed, the actual yield corresponds to 70% of the 
full phase space. A further increase towards equilibrium yield at higher RHIC- 
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200 energy range is expected. The remarkable feature of the RHIC situation is 
that much of strangeness enhancement is found in the (multistrange) baryon 
abundance. Given the large strangeness per baryon ratio, Fig.|l|, baryons and 
antibar5pLns produced from QGP at RHIC have been predicted to be mostly 
strange EJ. 

The high specific strangeness yield s/5 is a clear indicator for the extreme 
conditions reached in heavy ion collisions. An equally interesting observable 
is the occupancy of the strangeness phase space. In sudden hadronization, 
yHG/yQGP ~ 1^ the growth of volume is neghgible, T^^^ ~ T^^ , the tem- 
perature is maintained across the hadronization front, and the chemical occu- 
pancy factors in both states of matter accommodate the different magnitude 
of the particle phase space. In this case, the QGP strangeness when 'squeezed' 
into the smaller HG phase space results in 7^*^/7^*^^ — 3 , which is of the 
same magnitude as the unfrozen color degeneracy. This theoretical expec- 
tation is indeed observed both at top SPS and RHIC-130 run, as we have 
reported. 

For the top SPS energy range this interesting result could be ignored, 
since accidentally 7^*^/7^*^ — 1. Thus one can also model the hadronization 
at SPS energy in terms of an equilibrium hadronization model. The pion 
enhancement associated with the high entropy phase can be accommodated 
by use of two temperatures, one for the determination of absolute particle 
yields, and another for determination of the spectral shape. Such an approach 
has similar number of parameters, and comparable predictive power, the only 
inconsistency (with HBT) is the large volume required, as we have discussed. 
However, the condition, 7^^/7^'^ ~ 1, is not present at the RHIC energy 
range, where already at RHIC-130 the hadron phase space occupancy for 
strangeness is significantly larger than for light quarks, see table |[ It is the 
inclusion of the yields of multistrange antibaryons in the RHIC data analysis, 
which leads to this result. 

We see, at SPS and at RHIC, considerable convergence of the hadron pro- 
duction around properties of suddenly hadronizing entropy and strangeness 
rich QGP. The QGP phase strangeness occupancy rises from 7^*^^ ~ 0, 5 at 
SPS ysii? = 17.2 GeV, to 7^^^ ~ 0.7 at ^sSn = 130 GeV, and there 
is still space for a further strangeness yield rise at highest RHIC energy 
y^SNN = 200 GeV. In conclusion: the deconfined thermal QGP phase man- 
ifests itself through its gluon content, which generates in thermal collision 
processes a clear strangeness fingerprint of QGP. 
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